Candidate Pre-Mainsequence F Stars with Circumstellar Dust 
Identified Using Combined 2MASS and uvby Data 
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ABSTRACT 

We propose a method that uses near-infrared plus uvby photometry to identify po- 
tentially extensive circumstellar dusty environment about F and A stars. The method 
has been applied to a sample of ~ 900 metal rich reddened F stars with 2MASS and 
uvby data, suggesting the presence of circumstellar dust emitting in the near infrared 
for ~ 70 stars. The logTg — My diagram suggests that most, if not all, of them are 
likely pre-mainsequence (PMS). They seem to be consistent with being a continuation 
of the class of Herbig Ae/Be PMS stars into the spectral type F. Their number drops 
sharply downward of logTe ~ 3.84 (spectral types later than ~ F5), which may provide 
new clues to the PMS evolution of 1 — 2 Mq stars. We present a list of 21 most conspic- 
uous candidate stars with circumstellar dust. About half of them are associated with 
the extended star-forming region around p Oph. The brightest of these 21 stars, with 
V < 7.5, turn out to be IRAS sources, suggesting the presence of heated dust emitting 
in the far infrared. Also in this list, HD 81270 is reported as a very unusual star moving 
away from the Galactic plane at a projected speed of 70 kms~^. 

Subject headings: infrared: stars — planetary systems: protoplanetary disks — stars: circumstel- 
lar matter — stars: pre-main-sequence 



1. Introduction 

The detection of circumstellar dust disks 
provides us with an unprecedented opportu- 
nity to study the formation and evolution of 
exo-solar systems (see, e.g., Beckwith 1999) 
as well as the PMS stellar evolution (e.g., van 
den Ancker et al. 1997), given that many 
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of these systems are associated with PMS 
stars. The bulk of the known PMS stars are 
comprised by Herbig Ae and T Tauri stars. 
The presence of large amounts of dust about 
many of these stars has been inferred from 
IRAS observations (Aumann et al. 1984; Au- 
mann 1985; Gillett 1986; Backman & Gillett 
1987; Walker & Wolstencroft 1988; Backman 
& Paresce 1993). The IRAS survey still re- 
mains a major source to search for candi- 
date dust-disk systems (e.g., Mannings &: Bar- 
low 1998). Some 400 IRAS candidates ex- 
hibit an unusually large far-infrared (FIR) ex- 
cess (e.g., Stencel & Backman 1991, Man- 
nings & Barlow 1998), interpreted as being 
due to thermal emission of hot circumstellar 
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dust (30 < T < 1000 K). In contrast to IRAS 
systems whose large FIR excess results from 
a high-density concentration of warm dust, 
there are young stellar objects (YSO) that 
show relatively small FIR excess. They in- 
clude T Tauri stars exhibiting spectral energy 
distribution (SED) that declines longward of 
~ 2 /xm (YSOs with SED of Class II and Class 
III); these stars are interpreted as advanced 
evolutionary phases of T Tauri systems, with 
only a fraction of the original dust mass re- 
maining in the disk (e.g., Lada and Wilking 
1984; Andre & Montmerle 1994). Less mas- 
sive dust disks or envelopes may not be de- 
tectable in the far infrared. Disk clearing may 
have started for these stars, but it may still be 
possible to detect the disk remnant material 
in the near infrared through reflected light. 

In this paper we present a new criterion 
that uses a combined set of 2MASS and uvby 
data to identify F stars with circumstellar 
dust. The candidate stars found with this 
method appear to possess dusty disks or en- 
velopes radiating the reflected and/or thermal 
near-infrared (NIR), causing a star to look in 
the NIR brighter than expected for a reddened 
normal F star at the same visual magnitude 
V. Many of them are not IRAS sources and 
may therefore represent a dusty environment 
substantially different from that found among 
the IRAS sources. 

2. Data 

Starting with a sample described in Suchkov 
& McMaster (1999), which includes F stars 
with uvby colors from Hauck & Mermilliod 
(1998) and observed by the Hipparcos, we 
found ~ 5230 counterparts in the 2MASS 
Point Source Catalog. Some of the stars 
among those with high reddening are, in fact, 
of the spectral type A rather than F. The rea- 
son is that the selection for the original sample 
was made based on the {b — y) color, specifi- 
cally, 0.222 <h-y < 0.390. For unreddened 



stars, the latter criterion effectively selects 
stars later than spectral type A. However, if 
reddening is present, stars with intrinsically 
bluer colors, including A stars, can also be 
picked up by this criterion. The overall frac- 
tion of A stars in our sample is very small, and 
for the sake of simplicity, we'll be referring to 
the sample stars as F stars. The criterion also 
selects a small fraction of early type G stars. 

The uvby data were used to derive redden- 
ing Efj-y, effective temperature, and metal- 
licity as described in Suchkov & McMaster 
(1999) and Suchkov (2001), based on the algo- 
rithm in Moon (1985) as well as the Schuster 
& Nisscn (1989) mctallicity calibration. Ex- 
tinction was estimated as 



Av = 



i^.-.(3_354^0^040|E^^ 



0.74 



0.74 



(1) 



This equation follows from the results in 
Straizis (1977), Crawford (1975), and Sudzius 
(1974). Absolute magnitude was computed 
from the V magnitude corrected for extinc- 
tion, using the Hipparcos parallax; V magni- 
tudes were taken from the Hipparcos catalog. 

3. Circumstellar dust from intrinsic 
near-infrared excess 

The sample stars exhibit a tight temperature- 
induced correlation between the temperature- 
sensitive uvby color index (6 — y)o and V — K. 
In Figure 1, this correlation is shown for stars 
within 80 pc with small or no reddening, 
E})-y < 0.015. The scatter about the regres- 
sion line, ~ 0.05 after 3-sigma correction, is 
due to both photometric errors and real phys- 
ical effects, mainly related to absolute magni- 
tude, duplicity, and remaining reddening. It 
reduces to ay-K ~ 0.035 if the range for these 
parameters is narrowed down to 2 < My < 3, 
Eii-y < —0.005, and binaries are excluded. 
This should be close to the value resulting 
from photometric errors only. The linear re- 
gression shown in Figure 1 remains, however, 
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the same, within uncertainties, as for the con- 
strained sample. It was used to compute the 
deviation from the regression hne for ah sam- 
ple stars, 

5{V - K),_y = {V-K)-{V- K),_y. (2) 

A star having an excessive NIR flux with re- 
spect to the one expected for an unreddened 
standard star at the same effective tempera- 
ture would deviate upward from this line, 



d{V - K)b-y > 0. 



(3) 



Circumstellar dust can contribute to ob- 
served infrared excess in two quite different 
ways. First, this material would absorb stellar 
light stronger at shorter wavelengths, which 
adds up with the effect of interstellar extinc- 
tion and results in flux deficiency in the blue 
and visual with respect to the flux in the in- 
frared (the reddening effect). The part of the 
infrared excess due to this effect will be re- 
ferred to as apparent infrared excess. Second, 
a dust disk is a source of reflected and thermal 
radiation. Unlike circumstellar and interstel- 
lar extinction, this radiation actually adds to 
the obscured stellar flux in the infrared. The 
part of the near-infrared excess due to this 
extra flux will be referred to as intrinsic near- 
infrared excess. Isolating the latter part in 
the observed total infrared excess would thus 
amount to the detection of circumstellar dust. 

The two parts of the infrared excess can 
be separated as follows. The reddening effect 
in the total near-infrared excess as defined by 
equations (2) and (3) is clearly seen in the 
Eb-y — 6{V — K)i,^y diagram in Figure 2. In 
the diagram, the bulk of the reddened stars 
tightly envelope the line representing the re- 
lationship between 6{V — K)h-y and E^^y ex- 
pected from reddening only, 
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Fig. 1.- {V — K) versus (6 — y)o diagram 
for F stars within 80 pc. The thick line is 
the linear regression used to compute 5{y — 
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Fig. 2. — Near-infrared excess versus red- 
dening for metal-rich ([Fe/H] > 0) F stars. 
The most prominent candidate dust-disk stars 
are shown as filled circles, with IRAS sources 
among them indicated by overlaid asterisk. 
The solid line is the relationship between 
6{V — K)f,-y and E^-y expected from the red- 
dening effect only (see text for details). 



6{V - K)b-y = Av{Eb-y), 
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Fig. 3. — Comparison of stars with and 
without intrinsic near-infrared excess, A(y — 
K)i,^y > 0.07 (upper panel) and < -0.03 
(lower panel), respectively. Overlaid are 
isochrones for the Hyades chemical composi- 
tion (Chaboyer et al. 1999) labeled by age in 
Gyr. 



where Ay is defined by equation (1), and 
extinction in K is neglected in comparison 
with Ay- Suppose a star has a circumstel- 
lar dust disk reradiating stellar optical and 
UV light in the infrared as well as reflecting 
unabsorbed stellar infrared radiation. In the 
Ei,-y—d{V—K)b-y diagram, such a star would 
deviate upward from the line given by equa- 
tion (4), the amount of deviation representing 
a measure of intrinsic infrared excess. There- 
fore, we can use this deviation criterion 
identifying candidate dusty stars: 

A{V-K)h-y = S{V-K)h-y-Av{Eh-y) > 0. 

(5) 

We have applied this criterion to our sample 
of F stars. Given that large masses of dust are 
expected only in very young systems, it makes 
sense to start a search for dusty stars by ex- 
amining metal rich stars because high metal- 
licity typically implies a young age. There- 
fore, we have restricted our search sample to 
stars with [Fe/H] > 0. Also, since we expect a 
dusty system to be reddened, only stars with 
Eb^y > have been considered. These con- 
straints reduce the working sample to ~ 900 
stars. These stars are shown in Figure 2; 
those above the reddening line are predicted, 
in accordance with our criterion, to possess 
circumstellar dust. 

The most conspicuous candidate dusty 
stars are shown in Figures 2 and 3 as filled 
circles; they are also listed in Table 1. These 
21 stars have Eb_y > 0.01 and A{V-K)b-y > 
0.1. About half of them are spatially associ- 
ated with the giant molecular complex cen- 
tered at the site of ongoing star formation in 
p Oph. 

IRAS fluxes in Table 1, F(12 //), F(25 //), 
and F{60 /x), are from the IRAS Faint Source 
Catalog, Version 2. In total, we found seven 
out of the 21 selected candidate PMS dusty 
F stars to be IRAS sources (marked in Fig- 
ure 2 by overlaid asterisk) , which suggests the 
presence of heated circumstellar dust about 
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them. These 7 stars are brighter than V ~ 
7.5 mag (with a notable exception for HD 
144432, which is 8.188 in V), which may hint 
that the thermal emission from possible dust 
disks around fainter F stars is simply be- 
low the detection limit of the IRAS; however, 
an intrinsic difference between the IRAS and 
non-IRAS stars cannot be ruled out. We have 
also searched through the ISO archive (using 
SIMBAD interface) to see if any of the stars 
in Table 1 was among the ISO pointed obser- 
vations, but found none. 

To gain insight into the nature of the pre- 
dicted stars with circumstellar dust, we have 
examined the log Tg — My diagrams for more 
extreme, 2 a candidates, A(V — K)^p_y^^ > 
0.07, and compared it with the one for "non- 
candidates" defined by A{V — K)(^h-y)o < 
—0.03, selecting both groups from the sam- 
ple shown in Figure 2. Inspection of Fig- 
ure 3 reveals a large difference between the 
two diagrams. While the stars in the lower 
panel scatter across the entire temperature 
range displayed in the diagram, the stars in 
the upper panel exhibit a heavy concentra- 
tion at the hot end of the temperature range 
where only very young stars can reside. Thus, 
the stars predicted to possess dusty disks (up- 
per panel), are mostly very young whereas the 
non-candidates (lower panel) are dominated 
by a much older, on average, population. This 
is exactly what one would expect comparing 
stars with and without circumstellar dust. 

The prominent clustering of the stars in 
the upper panel in Figure 3 within a narrow 
temperature range of 3.86 < logTg < 3.83 
(see also Figure 4), which span about 3 mag 
in My, offers a further clue to the nature of 
these objects. The large spread of absolute 
magnitude they exhibit is obviously inconsis- 
tent with the effects of stellar evolution away 
from the main sequence; if it were, the dom- 
inant population in the log Tg — My diagram 
would have been at cooler rather than hotter 



temperatures, simply because the initial mass 
function increases toward lower masses. To 
put it differently, once we ask the question: 
"Where are stars in that diagram that must 
be numerous downward of log Tg ~ 3.83?" , we 
would immediately realize that most of the 
hot stars above the ZAMS in Figure 3 are 
not post-ZAMS stars but rather are at vari- 
ous phases of PMS evolution — which perfectly 
matches the inferred presence of dust around 
these stars. 

One may notice that also in the lower panel 
in Figure 3 and the respective middle panel in 
Figure 4, there are many hot stars separated 
from the bulk of cool stars by a poorly pop- 
ulated temperature interval around logTg ~ 
3.832. The same arguments as above make us 
believe that these are also mostly PMS stars. 
Additional argument is that their reddening is 
about twice as large as that of the cool stars, 
E{b-y) = 0.028 ±0.003 versus 0.012 ±0.001; 
this should be expected if a significant part of 
it, if not all, occurs in a dense circumstellar 
environment of a PMS star. The absence of 
excessive NIR emission apparently indicates 
that the dust is colder than in the PMS can- 
didates in the upper panel in Figures 3 and 4. 
Incidentally, these latter candidates are less 
reddened, E{b - y) = 0.021 ± 0.001, then the 
stars with presumably cold dust, which should 
be expected if warm dust becomes more dom- 
inant when much of the original dust content 
has been dispersed. 

Note that a large fraction of dusty PMS 
candidates inferred for the sample shown in 
Figure 3 is not at all typical for the general 
population of the field F stars. The temper- 
ature distribution of a sample that is more 
or less representative of a general population 
(the sample shown in Figure 1) is displayed 
in the lower panel of Figure 4. It is obvious 
that for this sample the fraction of stars in the 
temperature range dominated by PMS candi- 
dates is quite low and consistent with the ex- 
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Fig. 4. — Temperature distribution of the 
stars shown in Figure 3 (upper and middle 
panels) and Figure 1 (lower panel). Percent- 
age of stars to the right and left of the vertical 
line is indicated. 



pectations for an IMF declining toward higher 

stellar masses. 

The dusty PMS candidates are contiguous 
in the log Tg — My diagram with A stars that 
are known to contain a population of PMS 
emission-line, dusty stars usually referred to 
as a group as Herbig Ae/Be (HAeBe) stars. 
The properties of the candidates and the fact 
that they are adjacent to A stars suggest that 
they represent a continuation of the popula- 
tion of HAcBc stars into the spectral class F. 
The ct al. (1994) extended the criteria for 
membership of the class of HAeBe stars to 
include emission line objects of spectral types 
earlier than F8. However, the number of stars 
later than F2 in their catalog of HAeBe ob- 
jects is very small, only about half a dozen. 
The candidate PMS dusty F stars identified 
in this paper would represent an enormous in- 
crease in the number of F stars later than F2 
that appear to be akin to HAeBe stars. 

As discussed above, the stars in the upper 
panel of Figure 3 are concentrated within a 
narrow temperature range, 3.860 < logTg < 
3.832 (spectral range ^ F2 to F5), their num- 
ber plummeting downward of logTg ~ 3.84. 
This may have interesting implications perti- 
nent to the nature of the class of HAeBe ob- 
jects extended into the spectral range F, sug- 
gesting that such objects get extremely rare 
red ward of F5. This phenomenon is likely re- 
lated to fundamental properties of the PMS 
evolution of 1 — 2 Mq stars. So the detec- 
tion of the cut-off in the temperature (spec- 
tral type) distribution of HAeBe-|-F objects at 
logTe 3.84 (Sp ~ F5) should be expected to 
provide new insight into the physics of PMS 
stars. 

In summary, we have proposed a new crite- 
rion to identify F stars with warm circumstel- 
lar dust emitting in the near infrared. Appli- 
cation of this criterion to a sample of ~ 900 
metal rich, moderately reddened F stars re- 
sulted in detection of ~ 70 2-sigma candi- 
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dates. They appear to be PMS stars, rep- 
resenting a continuation of the class of Herbig 
Ae/Be stars into the spectral range F and ex- 
hibiting a cutoff redward of ~ F5. There is 
also evidence for a population of PMS stars 
with cold dust undetectable in the near in- 
frared. 
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Table 1 

Photometric properties of selected stars with intrinsic infrared excess 



Star 


Sp 


V 


J 


H 


K 


V-K 


F(12 m) 

(Jy) 


F(25 m) 

(Jy) 


F(60 At) 

(Jy) 


HD 18466 


A2/A3V 


6.286 






4.678 


1, 


.608 


0.555 


0.113 


0.112 


HD71766 


F2III 


5.995 






4.678 


1, 


.317 


0.495 


0.105 


0.113 


HD73378 


A9III 


8.000 


7.178 


6.997 


6.948 


1, 


.052 








HD81270* 


F6V 


8.950 


7.915 


7.636 


7.553 


1, 


.397 








HD88021 


F5 


6.678 


5.405 




4.915 


1, 


.763 


0.455 


0.231 


0.133 


HD115984 


F3V 


8.191 


7.405 


7.229 


7.121 


1, 


.070 








HD120850 


FOV 


7.362 


6.605 


6.449 


6.341 


1, 


.021 








HD121036 


F2V 


9.549 


8.609 


8.458 


8.370 


1, 


.179 








HD133411 


FOV 


7.911 


7.014 


6.888 


6.831 


1, 


.080 








HD144432'' 


A9/F0V 


8.188 


7.115 


6.541 


5.894 


2, 


.294 


7.62 


9.23 


5.68 


HD151070 


F5III 


6.871 


5.612 


5.297 


5.190 


1, 


.681 


0.389 


0.116 


0.131 


HD152500 


FOV 


8.202 


7.343 


7.227 


7.125 


1 


.077 








HD153777 


F2IV 


8.435 


7.547 


7.331 


7.224 


1, 


.211 








HD156441 


A9V 


9.100 


8.410 


8.245 


7.930 


1, 


.170 








HD158063 


F8 


7.383 


6.088 


5.745 


5.666 


1, 


.717 


0.219 


0.044 


0.092 


HD160410 


FOV 


8.210 


7.416 


7.109 


6.921 


1, 


.289 








HD 164200 


F3V 


8.410 


7.445 


7.306 


7.173 


1, 


.237 








HD172550 


FOV 


8.922 


8.116 


7.982 


7.862 


1, 


.060 








HD 195678 


A9V 


9.133 


8.297 


8.127 


8.077 


1, 


.056 








HD212500 


F4III 


6.950 


6.009 


5.863 


5.791 


1, 


.159 








HD222355 


F5 


7.548 


6.219 


5.907 


5.751 


1, 


.797 


0.204 


0.186 


0.173 



"Flagged in SIMBAD as high proper- motion. With tangential velocity of 101 km s~^, the star moves 
away from the Galactic plane at an extremely high velocity of 70 kms~^. The star is very unusual 
in other respects; it is extremely underluminous (below the Hyades ZAMS by ~ 1.2 mag), shows a 
signature of fast rotation and/or anomalously low surface gravity, and its metallicity exceeds the mean 
metallicity of the Hyades F stars by ~ 0.2 dex — all suggesting that something dramatic happened to 
this star recently. 

'^Flagged in SIMBAD as pre-mainsequence. Listed as a candidate disk-bearing system in Mannings 
& Barlow (1998) on the basis of its IRAS emission. 
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